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ABSTRACT 

A  series  of  flight  tests  were  conducted  on  the  PC  9/A  aircraft,  A23-045, 
at  the  Royal  Australian  Air  Force’s  Aircraft  Research  and  Development  Unit. 
System  identification  techniques  were  applied  to  the  data  obtained  firom  these 
flight  tests  to  determine  the  stability  and  control  derivatives  of  the  aircraft. 
The  longitudinal  results  for  the  aircraft  in  cruise  configuration  are  presented 
in  this  report  and  comparisons  are  made  with  empirical  and  ground  based 
estimates. 
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A  Correlation  between  Flight-determined  Longitudinal 
Derivatives  and  Ground-based  Data  for  the  Pilatus  PC  9/ A 
Training  Aircraft  in  Cruise  Configuration 


EXECUTIVE  SUMMARY 

Air  Operations  Division  (AOD)  has  developed,  or  acquired,  a  number  of  fixed-wing 
fiight  dynamic  models  and  is  also  responsible  for  providing  advice  to  the  Australian  De¬ 
fence  Organisation  (ADO)  on  flight  simulator  flight  dynamic  model  requirements.  The 
models  generally  make  use  of  extensive  static  and  dynamic  stability  and  control  derivative 
databases,  in  addition  to  engine  and  flight  control  models.  The  static  model  data  may  be 
obtained  firom  wind  tunnel  testing,  whilst  the  dynamic  data  is  traditionally  obtained  from 
flight  tests  using  system  identification  techniques. 

The  system  identification  techniques  used  to  estimate  aerodynamic  derivatives  of  con¬ 
ventional  aircraft  are  well  established.  The  major  requirement  of  these  techniques  is  high 
fidelity  measurements  of  manoeuvre  input  (i.e.  control  surface  deflections),  and  aircraft 
response  (i.e.  angular  rates  and  linear  accelerations),  as  well  as  air  data  including  airspeed, 
altitude,  angle-of-attack  and  angle-of-sideslip. 

Following  the  completion  of  the  AOD  PC  9/A  wind  tunnel  tests,  the  requirement  ex¬ 
isted  for  a  dedicated  flight  dynamic  modelling  flight  test  program  to  both  validate  the 
flight  dynamic  model  of  the  PC  9/ A  and  to  provide  dynamic  derivative  estimates  which 
were  unobtainable  in  the  AMRL  wind  tunnel.  An  instrumentation  suite,  including  an  air 
data  probe  for  the  direct  and  accurate  measmement  of  angle-of-attack,  angle-of-sideslip, 
temperature,  and  static  and  dynamic  pressure,  was  fitted  to  the  aircraft  and  a  test  ma¬ 
noeuvre  matrix  was  designed  specifically  for  the  purpose  of  gathering  flight  dynamic  data. 
This  flight  test  program  was  conducted  at  the  Aircraft  Research  and  Development  Unit 
(ARDU)  during  1998/99. 

This  report  details  the  analysis  of  the  longitudinal  manoeuvres  carried  out  with  the 
aircraft  in  the  cruise  configuration.  The  static  and  dynamic  derivatives  thus  obtained  are 
compared  with  wind  tunnel  estimates  as  well  as  a  number  of  empirical  estimates  obtained 
from  alternative  sources.  A  discussion  of  some  of  the  difficulties  encoimtered  during  the 
estimation  process  is  also  included.  The  data  obtained  from  these  tests  will  be  used  in 
the  development  of  the  PC  9/ A  flight  dynamic  model  for  the  purpose  of  enhancing  AOD 
support  for  the  PC  9/A  fleet,  including  possible  upgrad.es  to  the  part-task  trainer. 
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Notation 

Body  axes  linear  accelerations  (g) 

Reference  span  (10.124  m) 

Pitching  moment  coefficient 

r'  — 

'-'m  —  qSc 

Pitching  moment  coefficient  due  to  angle-of-attack 
C'ma  =  ^  (per  degree) 

Pitching  moment  coefficient  due  to  rate  of  change  of  angle-of-attack 
C'ma  =  (per  radian) 

Pitching  moment  coefficient  due  to  elevator  deflection 
Cms,  =  ^  (per  degree) 

Pitching  moment  coefficient  due  to  pitch  rate 
(per  radian) 

Normal  force  coefficient 
Cn  —  qs 

Normal  force  coefficient  due  to  angle-of-attack 
CiVa  =  ^  (per  degree) 

Normal  force  coefficient  due  to  rate  of  change  of  angle-of-attack 
C'^■d  =  (P®"^  radian) 

Normal  force  coefficient  due  to  elevator  deflection 
^  (P®^  degree) 

Normal  force  coefficient  due  to  pitch  rate 
(per  radian) 

Longitudinal  force  coefficient 


Vertical  force  coefficient 
Reference  chord  (1.65  m) 

Ratio  of  the  regression  mean  square  to  the  residual  mean  square 
Gravitational  acceleration  (9.81  m/s^) 

Roll  moment  of  inertia  (kg.m^) 

Pitch  moment  of  inertia  (kg.m^) 

Yaw  moment  of  inertia  (kg.m^) 

Cross  product  of  moment  of  inertia  (kg.m^) 
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Normal  force  (N) 
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Squared  multiple  correlation  coefficient 
r  Yaw  rate  (rad/s) 

S  Reference  area  (16.29  m^) 

V  True  velocity  (m/s) 
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Xc.g.  Longitudinal  c.g.  position  (m) 
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d  Time  derivative  of  angle-of-attack  (rad/s) 

P  Angle-of-sideslip  (°) 

6a  Aileron  deflection  (6a  =  6a^  —  6aji)  (°) 
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4>  Roll  angle  (°) 
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Subscripts 

b  Body  axes 

c.g.  Centre  of  gravity 
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L  Left  (port) 

R  Right  (starboad) 
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1  Introduction 

The  Pilatus  PC  9/ A  is  one  of  a  number  of  high  performance  turbo-prop  aircraft  currently 
operated  by  the  Royal  Australian  Air  Force  (RAAF).  Work  currently  tmderway  in  the  Air 
Operations  Division  (AOD)  of  the  Defence  Science  and  Technology  Organisation  (DSTO) 
includes  investigations  into  the  propeller  power  effects  of  such  aircraft,  and  the  PC  9/A 
was  considerd  a  suitable  platform  for  study.  A  six  degree-of-freedom  flight  dynamic  model 
of  the  PC  9/ A  has  been  developed  by  AOD  for  use  in  these  investigations,  as  well  as  for 
pilot-in-the-loop  simulations  in  the  Air  Operations  Simulation  Centre  and  incident  and 
accident  investigations. 

The  development  of  a  flight  dynamic  model  requires  information  on  the  static  and  dynamic 
stability  and  control  derivatives  of  the  aircraft,  as  well  as  flight  control  laws  and  physical 
properties.  The  static  data  for  the  PC  9/A  were  collected  during  both  power-off  and  power- 
on  testing  of  a  scaled  aircraft  model  in  the  Aeronautical  and  Maritime  Research  Laboratory 
(AMRL)  low-speed  wind  tunnel.  Additional  power-off  static  data  were  obtained  fi:om  a 
computational  fluid  dynamic  model  for  a  hmited  number  of  aircraft  configurations. 

Flight  test  data  for  the  PC  9/A  was  required  to  both  validate  the  flight  dynamic  model 
and  to  provide  dynamic  derivative  estimates  of  the  aircraft.  A  flight  test  program  was 
conducted  on  a  PC  9/A  aircraft,  serial  number  A23-045,  by  the  RAAF  Aircraft  Research 
and  Development  Unit  (ARDU)  between  November  1998  and  February  1999.  The  details 
of  this  flight  test  program  and  the  aircraft  instrumentation  are  reported  in  reference  [1]. 
The  estimation  of  the  stability  and  control  derivatives  from  the  flight  test  data  involved  the 
use  of  system  identification  techniques,  and  required  control  input,  aircraft  response  and 
flight  condition  data,  measured  using  a  high  fidelity  instrumentation  system.  Maximum 
likelihood  and  stepwise  regression  techniques  were  employed  for  the  system  identification. 

This  report  presents  the  longitudinal  stability  and  control  derivatives  of  the  PC  9/ A 
extracted  from  flight  test  data  for  the  cruise  configuration,  including  a  comparision  with 
AMRL  wind  tunnel  test  data  and  empirical  results.  Sections  2  and  3  present  details  of 
the  test  aircraft  and  instrumentation  system.  Section  4  presents  the  system  identification 
techniques,  while  section  5  discusses  the  results. 


2  PC  9/ A  Test  Aircraft 

2.1  Aircraft  Description 

The  flight  test  aircraft,  A23-045  is  a  Pilatus  PC  9/A  operated  by  ARDU.  The  PC  9/A  is 
a  single-engine,  metal-skinned,  low-wing,  tandem  two-seat  training  aircraft.  The  aircraft 
is  powered  by  a  Pratt  and  Whitney  PT6A-62  turbo-prop  engine  flat  rated  to  950  SHP  [2], 
which  drives  a  Hartzell  HC-D4N-2A  four-blade  variable  pitch  propeller.  The  aircraft  was 
instrumented  as  outlined  in  section  3. 
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2.2  Flight  Control  System 

The  aircraft  primary  flight  controls  consist  of  the  ailerons,  rudder  and  elevator.  The  control 
surfaces  axe  Tnfl.rma.11y  operated  from  a  conventional  dual  control  column  and  rudder  pedal 
arrangement.  The  stick  and  rudder  pedals  are  connected  to  the  control  surfaces  through 
a  system  of  control  rods,  bellcranks,  cables  and  levers.  Trimming  control  is  provided  on 
aU  three  axes. 


2.3  Weight,  Centre-of-Gravity  and  Mass  Moments-of-Inertia 

The  test  aircraft,  A23-045,  was  weighed  by  ARDU  prior  to  the  flight  tests  and  had  a  basic 
mflfls  of  1784.5  kg  and  a  longitudinal  centre-of-gravity  position  of  26.25%MAC  [3].  Dining 
the  flight  test  program,  the  aircraft  weight,  centre-of-gravity  and  mass  moments-of-inertia 
varied  with  fuel  usage.  These  parameters  were  calculated  for  each  test  manoeuvre  using 
the  equations  in  Appendix  A. 


3  Instrumentation 

Aircraft  A23-045  was  fitted  with  an  instrumentation  system  designed  specifically  for  the 
gathering  of  flight  dynamic  data.  A  summary  of  the  instrumentation  used  in  the  flight 
test  program  is  included  below  and  a  more  comprehensive  description  of  the  design  re¬ 
quirements  and  calibration  is  included  in  references  [1]  and  [4].  The  ARDU  General  Data 
Acquisition  System  (GDAS)  was  fitted  in  place  of  the  rear  ejection  seat.  The  data  were 
encoded  by  a  16-bit  pulse  code  modulation  system  and  were  recorded  onboard  the  air¬ 
craft  using  a  MARS-2000  14-track  tape.  Real-time  flight  test  monitoring  was  provided  by 
telemetry  data  transmitted  to  the  ARDU  Primary  Analysis  Processor  (PAP)  hut. 

The  angular  rates  ip,q,r)  and  linear  accelerations  {axiCiY,CLz)  were  measured  using  the 
ARDU  KAISG 1134-1  Motion  Platform.  This  consisted  of  three  Smith  Industries  950  RGS 
flTign1fl.r  rate  gyros  and  three  SunStrand  QA1400  servo  accelerometers.  The  aircraft  roll, 
pitch  and  yaw  attitude  angles  (0,  6,  tj})  were  obtained  by  tapping  output  from  the  existing 
LISA  2000A  Artificial  Horizon  Reference  System  (AHRS). 

Static  outside  air  temperatiure  (OAT),  indicated  airspeed  (IAS),  angle-of-attack  (a)  and 
angle-of-sideslip  (j0)  were  obtained  from  the  Rosemount  Model  92AN  flight  test  air  data 
boom,  mounted  on  the  outboard  hardpoint  on  the  starboard  wing  (see  figure  1).  Aileron 
and  elevator  deflections  were  measured  using  Space  Age  Control  Inc.  series  160  cable 
position  transducers.  Rudder  deflection  was  measured  using  a  type  26V-11CX4C  position 
transducer.  All  sensors  were  calibrated  prior  to  commencement  of  the  flight  test  program. 

The  aircraft  was  also  instrumented  to  measure  engine  torque,  propeller  speed,  inlet  turbine 
temperature  (ITT),  gas  generator  speed  (NG),  fuel  flow  and  fuel  quantity.  The  torque, 
propeller  speed,  altitude  and  true  airspeed,  when  used  in  conjunction  with  a  performance 
map  for  the  Hartzell  HC-D4N-2A  propeller  [5],  allowed  the  calculation  of  engine  thrust. 
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4  Methods  of  Analysis 

During  the  flight  test  program,  doublet  and  3-2- 1-1  manoeuvres  were  performed  about 
steady  flight  conditions  at  airspeeds  between  90  and  200  KIAS.  Reference  [1]  describes 
these  manoeuvres  in  detail.  The  34  longitudinal  doublet  and  140  longitudinal  3-2-1-1  ma¬ 
noeuvres  were  analysed  using  stepwise  regression  [6]  and  maximum  likelihood  estimation 
techniques  [7]  to  determine  the  longitudinal  stability  and  control  derivatives  of  the  aircraft. 

Aircraft  derivatives  estimated  from  AMRL  power-off  wind  tunnel  tests  [8]  and  empirical 
data  estimated  by  Pilatus  [9]  using  Digital  Datcom  [10]  [11]  were  used  to  provide  compar¬ 
isons  and  in  the  case  of  the  maximum  likelihood  analysis,  to  provide  a  priori  estimates  to 
increase  the  rate  of  convergence  of  the  system  identification. 

A  right  handed  orthogonal  axes  system  was  adopted  for  the  analysis  of  the  flight  test  data. 
Positive  control  surface  deflections  were  defined  as  elevator  trailing  edge  down,  rudder 
trailing  edge  to  port,  and  starboard  aileron  trailing  edge  down,  port  aileron  trailing  edge 
up. 


4.1  Stepwise  Regression 

Stepwise  regression  is  an  unbiased  least  squares  estimator  in  which  new  independent  vari¬ 
ables  are  inserted  into  a  model,  one  at  a  time,  imtil  the  regression  equation  is  deemed 
acceptable.  The  appropriateness  of  the  model  can  be  determined  by  examining  a  num¬ 
ber  of  quantities  including  the  squared  multiple  correlation  coeflicient  and  the  F  statistic. 
The  squared  multiple  correlation  coefficient,  gives  a  measure  of  the  importance  of 
each  variable  as  it  is  inserted  into  the  equation  [6];  however,  the  improvement  in  due 
to  the  addition  of  new  terms  must  have  some  real  significance  besides  simply  reflecting  the 
inclusion  of  more  terms.  This  can  be  determined  by  monitoring  the  F  statistic,  the  ratio  of 
the  regression  mean  square  to  the  residual  mean  square.  The  inclusion  of  any  significant 
terms  is  generally  accompanied  by  an  increase  in  the  F  statistic  and  the  best  fit  with  the 
least  niunber  of  parameters  may  be  obtained  by  maximising  F.  Any  variable  which  does 
not  make  a  significant  contribution  is  removed  from  the  model,  with  the  selection  process 
continuing  until  no  new  variables  remain  to  be  inserted  into  the  equation.  Whilst  stepwise 
regression  gives  estimates  of  the  derivatives  included  in  the  regression  equation,  it  also 
permits  a  suitable  structure  for  the  model  to  be  determined.  In  addition,  it  provides  an 
independent  check  on  the  data  estimated  using  the  maximum  likelihood  technique. 

The  stepwise  regression  technique  was  applied  using  code  available  in  the  MATLAB® 
Statistics  Toolbox  [12].  The  following  model  equations  were  identified  during  the  analysis. 


Cn  =  Cno  +  C'jVaQ!  +  Cns^  Se 


(1) 


Cm  —  Cmo  "b  Cma  ^  ' 


(2) 
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4.1.1  Error  Band 

Included  in  the  stepwise  regression  analysis  is  the  calculation  of  an  error  band  on  the 
estimated  derivatives.  For  a  confidence  interval  of  95%,  this  error  band  is  approximately 
equal  to  two  standard  deviations.  Figures  2  to  11  show  the  stepwise  regression  derivative 
estimates,  including  the  calciQated  error  band. 

4.2  Maximum  Likelihood 

The  mavirnmn  likelihood  estimation  technique  was  applied  using  the  computer  program, 
pEst,  developed  at  NASA  Dryden  Flight  Research  Center  [13].  pEst  is  an  interactive  pa¬ 
rameter  estimation  program  which  solves  a  vector  set  of  time- varying,  ordinary  differential 
equations  of  motion. 

The  longitudinal  derivatives  identified  using  the  maximum  likelihood  technique  are  given 
in  table  1.  Three  aerodjmamic  derivatives  used  in  the  longitudinal  flight  dynamic  model 
were  not  identified  due  to  both  their  small  contribution  and  difficulty  of  estimation. 


Normal  Force 

Pitching  Moment 

Aerodynamic 

Cnc 

Cma 

Cm. 

Control 

CNh. 

Cms. 

Table  1:  Longitudinal  derivatives  estimated. 

The  derivative  Cma  is  difficult  to  estimate  from  standard  flight  test  manouevres  and  is 
usually  considered  as  an  additional  component  of  the  pitch  damping  derivative  Cmq  during 
parameter  identification  [14].  Experience  gained  during  the  analysis  of  F-lllC  flight  test 
data  [15]  has  shown  that  it  is  preferable  to  fix  Cma  a*  it’s  a  priori  value,  chosen  here  as 
-7.823  (per  radian),  and  allow  Cm,  to  be  estimated. 

Difficulties  also  arise  when  trying  to  estimate  both  the  normal  force  due  to  pitch  rate 
derivative,  Cjv,,  and  the  normal  force  due  to  the  time  rate-of-change  of  angle-of-attack 
derivative,  Cjv^  •  In  the  case  of  the  PC  9/ A,  analysis  showed  that  the  relative  contributions 
of  Cjv,  and  Cn^  to  the  total  aircraft  normal  force  were  negligible,  and  these  derivatives 
could  be  fixed  at  their  a  priori  values,  chosen  as  7.96  (per  radian)  and  3.061  (psr  radian), 
respectively,  without  any  loss  to  the  accuracy  of  the  estimation  of  the  other  normal  force 
derivatives,  and  Cjyjg  (refer  appendix  B). 

The  a  priori  values  of  Cm^t  C^q  nnd  Cjv^  were  obtained  from  empirical  data  in  refer¬ 
ence  [9].  The  resulting  total  force  and  moment  coefficient  equations  used  in  the  state 
and  response  equations  for  the  maximum  likelihood  analysis  of  the  PC  9/A  are  given  in 
equations  3  and  4.  The  differences  between  equations  1  and  2,  and  3  and  4  are  important, 
and  arise  from  the  inclusion  of  the  a  priori  values  in  the  maximum  likelihood  technique. 
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nr  0[C 

Cm  =  Cmo  +  Cm^a  +  Cm,^  +  Cms,  4  "  7-823—  (4) 

4.2.1  Cramer-Rao  Bounds 

For  the  estimated  parameters,  pEst  calculates  a  measure  of  the  estimation  certainty  known 
as  the  Cramer-Rao  bound.  A  detailed  interpretation  of  this  quantity  is  given  in  [16].  The 
Cramer-Rao  bounds  are  shown  for  each  derivative  estimated  by  pEst  in  figures  2  to  11. 
The  Cramer-Rao  bounds  have  been  factored  in  accordance  with  the  procedures  described 
in  [16]  to  account  for  the  presence  of  band-limited  noise. 


5  Results  and  Discussion 

Longitudinal  derivatives  estimated  via  maximum  likelihood  parameter  estimation,  stepwise 
regression,  empirical  methods  and  power-off  wind  tunnel  experiments  are  summarised  in 
figures  2  to  11. 


5.1  Angle-of- Attack  Derivatives 

Figures  2  to  6  show  the  angle-of-attack  derivatives  and  Cma  plotted  against  oc.  Flight 
test  estimates,  derived  from  3-2-1-1  manoeuvres,  of  the  normal  force  due  to  angle-of-attack 
derivative,  C'^^ ,  display  very  good  agreement  with  the  wind  tunnel  estimate  from  both 
system  identification  techniques.  The  results  exhibit  a  similar  variation  with  angle-of- 
attack,  with  CNa  decreasing  as  a  increases,  as  that  predicted  by  wind  tunnel  tests.  For 
the  doublet  manoeuvres,  the  maximum  likelihood  estimates  fall  between  the  empirical  and 
wind  tmmel  estimates,  whilst  the  stepwise  regression  estimates  still  agree  well  with  the 
wind  tunnel  data. 

The  fiight  test  estimates  of  the  pitch  stiffness  derivative,  Cma,  generally  smaller  in 
magnitude  than  the  wind  tunnel  and  empirical  estimates  which  are  both  power-off.  This 
difference  is  due  to  the  imbalance  between  two  competing  effects  acting  on  the  portion 
of  the  horizontal  tail  within  the  propeller  slipstream,  these  being  the  increased  dynamic 
pressure  and  the  reduced  angle  of  incidence.  The  reduction  in  angle  of  incidence  dominates, 
resulting  in  an  overall  reduction  in  the  pitch  stiffness  of  the  aircraft,  as  seen  in  figures  2  to  6. 
This  phenomenon  is  discussed  further  in  section  5.4.  The  fiight  test  results  exhibit  the 
same  variation  with  angle-of-attack  as  predicted  by  the  wind  tunnel. 

5.2  Pitch  Rate  Derivatives 


Figures  2  to  6  show  the  pitch  rate  derivative  Cm,  plotted  against  a.  As  discussed  earlier, 
Cm^  was  fixed  at  it’s  a  priori  value  to  aid  the  estimation  of  the  pitch  damping  deriva¬ 
tive  Cmq-  This  method  has  been  used  successfully  during  other  parameter  estimation 
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exercises  [15],  and  the  results  shown  in  figures  2  to  6  confirm  this.  In  general,  the  fiight 
test  estimates  of  Cm,  from  maximum  likelihood  and  stepwise  regression  agree  well  with 
empirical  data.  The  AMRL  low-speed  wind  tunnel  has  no  facility  to  measure  dynamic 
derivatives.  Thus  only  empirically  derived  estimates  are  available  for  comparison.  The 
flight  test  estimates  exhibit  a  degree  of  scatter  at  angles-of-attack  higher  than  9°. 

Data  analysis  carried  out  previously  on  the  PC  9/A  and  F-lllC  [15]  has  consistently 
highlighted  the  difficulty  of  obtaining  reasonable  estimates  for  the  normal  force  due  to 
pitch  rate  derivative,  Cjv,.  In  the  past,  the  derivative  Cjv,  was  usually  fixed  at  some  a 
priori  value  whilst  the  remaining  derivatives  were  estimated.  The  opportunity  was  taken 
with  the  PC  9/A  data  to  further  investigate  the  effect  that  fixing  Cjv,  would  have  on  the 
other  derivatives.  This  investigation  found  that  the  contributions  of  both  Cjv,  and  to 
the  total  normal  force  coefficient,  Cn,  are  negligible,  hence  the  effect  of  any  errors  present 
in  the  a  priori  estimate  of  should  be  minimal.  A  detailed  discussion  of  this  analysis 
is  given  in  appendix  B. 


5.3  Control  Derivatives 

Figures  7  to  11  show  the  control  derivatives  and  plotted  against  a.  The 
Tnayimiim  likelihood  estimates  of  agree  well  with  the  wind  tunnel  estimate,  whilst 
stepwise  regression  consistently  imderestimates  the  magnitude  of  Cms^  in  comparison  with 
the  wind  tunnel  estimate.  The  fiight  test  results  also  exhibit  a  reasonable  degree  of  scatter 
above  9°  angle-of-attack.  Overall,  however,  the  low  Cramer-Rao  bounds  would  indicate  a 
high  level  of  confidence  in  the  maximum  likelihood  estimates  of  Cms^  • 

Estimates  of  show  a  high  degree  of  scatter  at  angles-of-attack  above  9°  and  fall  short 
of  both  the  wind  tunnel  and  empirical  estimates.  In  the  case  of  the  stepwise  regresion 
results,  negative  values  of  were  observed.  As  discussed  previously,  a  priori  values  are 
included  in  the  maximum  likelihood  technique,  whereas  the  stepwise  regression  technique 
does  not  consider  these  values.  In  order  to  assess  the  effect  of  this  difference,  a  number 
of  manoeuvres  were  analysed  using  the  maximum  likelihood  technique  where  the  a  priori 
values  were  set  to  zero.  This  analysis  showed  the  effect  of  the  a  priori  values  on  the 
estimated  derivatives  to  be  small.  However,  a  change  in  sign  of  was  observed,  similar 
to  that  seen  in  the  stepwise  regression  results.  It  is  therefore  proposed  that  the  effect  of 
setting  both  ^Jid  Cjv,  to  zero  is  to  reduce  the  rate  of  onset  of  normal  force  during  a 
manoeuvre  and  therefore,  the  modelling  techniques  decrease  the  value  of  ,  to  below 
zero,  to  compensate.  This  ensures  the  overall  match  of  the  data  is  good  but  results  in 
some  loss  of  model  resolution  during  the  inital  apphcation  of  elevator. 


5.4  Variation  of  Longitudinal  Derivatives  with  Rate  of  Climb 
and  Angle-of- Attack 

Figure  12  shows  the  variation  of  the  longitudinal  derivatives  with  angle-of-attack  and 
nominal  rate  of  climb  for  a  selection  of  3-2- 1-1  manoeuvres.  A  linear  curve  fit  was  applied 
to  each  rate  of  climb  data  set,  to  highlight  the  trends  shown  by  the  data.  The  figure  shows 
that  the  influence  of  rate  of  climb  on  the  longitudinal  derivatives  increases  with  increasing 
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angle-of-attack,  with  the  magnitude  of  Cma  decreasing  (an  effective  forward  movement 
of  the  aircraft  neutral  point),  while  the  magnitudes  of  Cmq  and  Cmsq  increase  (due  to 
the  increased  dynamic  pressure  in  the  propeller  slipstream).  These  trends  result  from  the 
increase  in  aircraft  thrust  coefficient  with  rate  of  climb  at  a  given  angle-of-attack.  They  are 
in  agreement  with  previous  theoretical  and  experimental  studies  detailed  in  references  [17] 
and  [18].  For  the  manoeuvres  in  figure  12,  thrust  coefficient  varies  fi:om  approximately 
0.02  to  approximately  0.09. 


5.5  Manoeuvre  Type  and  Altitude  Dependency 


Tables  2  and  3  show  the  mean  values  across  the  linear  angle-of-attack  range  (0°  -  6°) 
of  the  aerodynamic  derivatives  estimated  by  both  the  maximum  likelihood  and  stepwise 
regression  techniques,  for  each  manoeuvre  type  and  altitude  band. 


C'jVa 

Cns. 

Cm. 

Cma 

Cms. 

3211  5000  ft 

8.94E-02 

5.50E-03 

-7.70E-03 

-1.44E-h01 

-2.15E-02 

3211  10  000  ft 

8.82E-02 

5.40E-03 

-7.20E-03 

-1.41E-f01 

-2.12E-02 

3211  15  000  ft 

8.77E-02 

5.00E-03 

-7.00E-03 

-1.47E-H01 

-2.12E-02 

Doublet  5000  ft 

9.18E-02 

2.96E-03 

-7.74E-03 

-1.38E-t-01 

-2.10E-02 

Doublet  15  000  ft 

8.99E-02 

3.60E-03 

-7.18E-03 

-1.59E-1-01 

-2.18E-02 

Table  2:  Summary  of  derivatives  estimated  using  maximum  likelihood. 


Cn. 

Cm. 

Cm„ 

Cms. 

3211  5000  ft 

7.78E-02 

-3.00E-03 

-5.50E-03 

-1.93E-H01 

-1.98E-02 

3211  10  000  ft 

8.45E-02 

-4.30E-03 

-5.50E-03 

-1.39E+01 

-1.67E-02 

3211  15  000  ft 

7.93E-02 

-4.40E-03 

-4.80E-03 

-1.54E+01 

-1.61E-02 

Doublet  5000  ft 

8.14E-02 

-4.90E-03 

-6.43E-03 

-1.45E+01 

-1.61E-02 

Doublet  15  000  ft 

8.27E-02 

-4.44E-03 

-6.21E-03 

-1.64E-I-01 

-1.73E-02 

Table  3:  Summary  of  derivatives  estimated  by  stepwise  regression. 

The  flight  test  estimates  for  the  3-2- 1-1  and  doublet  manoeuvres  do  not  show  any  signifi¬ 
cant  dependency  on  manoeuvre  type.  In  general,  the  results  show  increasing  uncertainty 
bounds  and  scatter  with  increasing  angle-of-attack.  This  trend  was  more  evident  in  the 
3-2-1- 1  manoeuvres,  which  were  performed  over  a  larger  angle-of-attack  range.  This  is  not 
suprising,  as  at  trim  angles-of-attack  above  approximately  9°  the  aircraft  was  operating  in 
the  non-linear  region  of  the  flight  envelope.  Furthermore,  the  test  pilot  reported  difficulty 
in  maintaining  trim  and  performing  the  manoeuvres  at  these  conditions. 

As  all  the  derivatives  are  non-dimensionalised  by  dynamic  pressmre,  they  were  expected 
to  be  insensitive  to  the  altitude  at  which  the  manoeuvre  was  executed.  Figures  2  to  11 
confirm  this. 
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6  Conclusions 

Derivatives  describing  the  longitudinal  stability  and  control  characteristics  of  the  PC  9/A 
in  the  cruise  configuration  were  determined  firom  flight  test  data  measurements  using 
TnaYirmiTn  likelihood  and  stepwise  regression  techniques.  The  results  were  compared  with 
power-off  wind  tunnel  and  empirical  values. 

Estimates  of  the  longitudinal  stability  derivatives,  Cj^a »  ^ma  Cm, ,  showed  reasonable 
agreement  with  the  corresponding  wind  tunnel  values  and  empirical  estimates.  The  flight 
test  and  wind  tunnel  results  also  exhibit  similar  trends,  when  plotted  against  angle-of- 
attack.  Good  agreement  with  wind  tunnel  and  empirical  estimates  was  observed  for  the 
longitudinal  control  derivate,  Cmsg  •  The  results  for  showed  a  high  degree  of  scatter 
and  were  affected  by  the  exclusion  of  a  priori  values  in  the  stepwise  regression  analysis. 
Any  discrepancies  in  this  derivative  should  not  affect  the  aircraft  model  due  to  its  low 
significance  in  the  total  Cn  equation. 

In  general,  the  flight  test  derivative  estimates  did  not  show  any  significant  dependence 
on  altitude  or  manoeuvre  type.  However,  they  showed  larger  scatter  and  possessed  larger 
uncertainty  botmds  at  angles-of-attack  above  9°,  where  the  aircraft  was  operating  in  the 
non-linear  region  of  the  flight  envelope.  The  longitudinal  derivatives  Cm^ ,  Cm,  and  Cms^ 
exhibited  trends  when  plotted  as  a  function  of  climb  rate  that  were  consistent  with  those 
predicted  by  earlier  theoretical  and  experimental  studies. 
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Figure  2:  PC  9/A  angle-of-attaek  and  pitch  rate  derivatives,  3-2-1-1  manoeuvres  at  5000  ft. 
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Figure  4-  PC  9/A  angle-of-attack  and  pitch  rate  derivatives,  3-2-1-1  manoeuvres  at  15  000  ft. 
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Figure  6:  PC  9/A  angle- of- attack  and  pitch  rate  derivatives,  doublet  manoeuvres  at  15  000  ft. 
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Figure  7;  PC  9/A  control  derivatives,  3-2-1-1  manoeuvres  at  5000  ft. 
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Figure  8:  PC  9/A  control  derivatives,  3-2-1- 1  manoeuvres  at  10  000  ft. 
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Figure  9:  PC  9/A  control  derivatives,  3-2-1-1  manoeuvres  at  15  000  ft. 
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Figure  11:  PC  9/ A  control  derivatives,  doublet  manoeuvres  at  15  000  ft. 


21 


(per  degree)  0^  (per  radian)  Cmg  (per  degree) 


DSTO-TR-0937 


Appendix  A:  Weight,  Centre-of- Gravity  and 
Mass  Moments-of-Inertia 


The  aircraft  mass,  centre-of-gravity  and  mass  moments-of-inertia  were  determined  for  each 
test  manoeuvre  based  on  the  aircraft  fuel  mass.  The  basic  aircraft,  pilot  and  fuel  masses 
and  moment  arms  from  references  [2]  and  [3]  are  given  in  table  Al.  The  moment  arms  are 
measured  relative  to  the  aircraft  datum,  located  3m  forward  of  the  engine  firewall,  and 
2m  below  the  forward  fuselage  reference  plane. 


Mass  (kg) 

Moment  Arm  (m) 

Basic  Aircraft 

1784.5 

4.306 

Pilot 

81.6 

4.061 

Fuel 

Variable 

4.178 

Table  Al:  Flight  test  aircraft  mass  distribution  [2,  3]. 

The  mass  of  the  aircraft  is  calculated  as  the  sum  of  the  basic  aircraft,  pilot  and  remaining 
fuel  masses,  by  the  following  equation. 


M  =  1866.1  +Mf  (kg)  (Al) 

The  longitudinal  centre-of-gravity  position  of  the  flight  test  aircraft  was  determined  from 
the  following  equation. 


^c.g. 


_  8015.5+4.178AfF 
~  1866.1-(-Mf 


(M) 


(A2) 


During  the  flight  test  program,  the  aircraft  longitudinal  centre-of-gravity  position  varied 
between  4.27  m  (24.2%MAC)  and  4,30  m  (25.6%MAC).  The  lateral  and  vertical  centre- 
of-gravity  positions  were  assumed  to  be  invariant  with  fuel  usage,  and  were  fixed  at  values 
of  0.024m  and  -2.2  m,  respectively,  relative  to  the  aircraft  datum. 

The  mass  moments-of-inertia  of  the  aircraft  are  given  by  the  following  equations. 


Ixx 

=  2505.9 -b  6.177  (M- 

-  1866.1) 

(kg.m^) 

(A3) 

lyy 

=  6622.2  +  0.033  (M  - 

-  1866.1) 

(kg.m2) 

(A4) 

Izz 

=  8467.1  4- 6.188  (M- 

-  1866.1) 

(kg.m^) 

(A5) 
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IXY 

=  49.0 +  0.0057  (M- 

-  1866.1) 

(kg-m^) 

(A6) 

Ixz 

=  196.9  +  0.0041  (M  ■ 

-  1866.1) 

(kg.m^) 

(A7) 

Iyz 

=  3.0  + 0.0007  (M- 

1866.1) 

(kg.m2) 

(A8) 
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Appendix  B:  C^q  Contribution  to  Total  Cn 

Equation 

Data  analysis  carried  out  previously  on  the  PC  9/ A  and  other  aircraft,  has  consistently 
highlighted  the  difficulty  of  obtaining  reasonable  estimates  for  the  normal  force  due  to 
pitch  rate  derivative,  Cn^-  In  the  past,  the  derivative  Civ,  was  usually  fixed  at  some 
a  priori  value  whilst  the  remaining  derivatives  were  estimated.  The  opportunity  was 
taken  to  further  investigate  the  effect  that  fixing  would  have  on  the  other  derivatives. 
This  investigation  was  done  in  two  parts.  First,  the  total  contribution  of  the  pitch  rate 
derivative  to  the  total  normal  force  coefficient  was  determined.  Then,  the  effect  of  changes 
in  Cat,  on  the  Tna.yimiiTn  likelihood  estimates  of  the  remaining  derivatives  was  examined. 

The  normal  force  coefficient  may  be  represented  by  the  following  equation: 


Cn  =  Cno  +  Cnc  oi  +  CNg  ^  +  Cns^  (Bl) 

The  contribution  of  the  pitch  rate,  g,  on  the  normal  force  coefficient  of  the  PC  9/A  was 
determined  by  comparing  the  relative  sizes  of  each  component  in  the  above  equation  for  a 
sample  case.  Those  results  are  shown  in  figure  B1  and  it  can  be  seen  that  changes  in  the 
normal  force  coefficient  are  mainly  due  to  changes  in  the  aircraft  angle-of-attack  and  that 
the  contribution  of  the  pitch  rate  is  relatively  small. 

The  effect  of  error  in  the  a  priori  estimate  of  Cjv,  was  investigated  by  changing  the  value 
of  CNg  and  examining  the  changes  in  the  other  derivatives.  Results  from  the  sample  case 
are  shown  in  table  Bl  and  they  show  that  changes  in  CNg  appear  to  have  little  effect  on 
most  of  the  derivatives  with  the  exception  of  Cns^-  This  derivative  shows  considerable 
variation  with  CNg,  which  was  to  be  expected.  Figure  Bl  shows,  however,  that  like  the 
CNg  derivative,  Cn^^  has  only  a  small  contribution  to  the  total  normal  force  coefficient. 
Hence  the  effect  of  any  errors  present  in  the  a  priori  estimate  of  Cjv,  should  be  minimal. 


CNg 

-16.0 

-8.00 

0.00 

7.96 

8.00 

16.00 

Cno 

CNcc 

Cm, 

8.44E-02 

8.56E-02 

-l.llE-02 

9.56E-02 

8.39E-02 

-8.39E-03 

1.07E-01 

8.24E-02 

-8.24E-03 

1.18E-01 

8.08E-02 

-4.31E-03 

1.18E-01 

8.08E-02 

-4.50E-03 

. . 1 

1.29E-01 

8.79E-02 

-4.32E-03 

Table  Bl:  Effect  of  CNg  variation  on  Cn  equation  components. 
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Figure  Bl:  Relative  magnitudes  ofCf/  equation  components. 
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